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We perform a number of improvements to the previous AKK extraction of fragmentation functions 
for TT*, -ftT^, p/p, Kg and A/A particles at next-to-leading order. Inclusive hadron production 
measurements from pp{p) reactions at BRAIfMS, CDF, PHENIX and STAR are added to the data 
sample. We use the charge-sign asymmetry of the produced hadrons in pp reactions to constrain 
the valence quark fragmentations. Data from e'^e~ reactions in regions of smaller x and lower y/s 
are added. Hadron mass effects are treated for all observables and, for each particle, the hadron 
mass used for the description of the e'^e~ reaction is fitted. The baryons' fitted masses are found to 
be only around 1% above their true masses, while the values of the mesons' fitted masses have the 
correct order of magnitude. Large x resummation is applied in the coefficient functions of the e"'"e~ 
reactions, and also in the evolution of the fragmentation functions, which in most cases results in a 
significant reduction of the minimized . To further exploit the data, all published normalization 
errors are incorporated via a correlation matrix. 

I. INTRODUCTION 

In perturbative QCD, fragmentation functions (FFs) D'^{x, Al'j), which can be interpreted as the probability for 
a parton i at the factorization scale Mf to fragment to a hadron h carrying away a fraction x of its momentum, 
are a necessary ingredient in the calculation of single hadron inclusive production in any reaction. Interest in FFs is 
widespread, to be found for example in the study of the proposed hot quark gluon plasma (QGP) of the early universe 
currently being sought in heavy ion collisions, in investigating the origin of proton spin, and in tests of QCD such as 
theoretical calculations for recent measurements of inclusive production in pp collisions at RHIC. 

From the factorization theorem, the leading twist component of any single hadron inclusive production measurement 
can be expressed as the convolution of FFs, being the universal soft parts containing the final state, with the equivalent 
productions of real partons, which are perturbatively calculable, up to possible parton distribution functions (PDFs) 
to account for any hadrons in the initial state. Thus, by using these data to constrain the FFs, predictions for future 
measurements can be made from current data. An exception to this possibility occurs when some new measurement 
depends on any regions of the FFs' function space that has not yet been constrained by experiment. For this reason, 
a failure to describe some data does not imply irrefutably that there are relevant physics effects which have been 
neglected in calculations. In particular, the apparent inconsistencies within this framework of universality occurring 
between charge-sign unidentified hadron production in e+e~, for which the contribution from gluon fragmentation is 
much less than from quark, and in pp(j)) reactions may be attributed to the large experimental uncertainties on the 
gluon fragmentation determined from e~^e~ reaction data only, rather than to neglected effects in the description of 
both reactions such as higher twist, heavy quark masses, resummation at large and small x, and higher order terms in 
the perturbative approximation, or to the less well understood hot QGP invoking parton energy loss. A combined fit 
of FFs to data from e~^e~ and pp{p) reactions would prove the optimum method of verifying consistency between the 
two types of reactions, since a successful fit to both types of data would imply that these apparent inconsistencies in 
fact lie within the experimental and theoretical uncertainties and so are not inconsistencies at all. Success is expected 
for identified particles in general, since good agreement is found [ij for the theoretical calculation for pp ^ tt + X 
data, where tt = 7r° 0] and [s], using FFs for n"^ constrained by data for e+e^ + X processes, and data 

for the production of tt^ is generally more accurate and plentiful than for the production of other particles due to 
the high abundance of tt^ in the particle sample. In other words, the current theoretical state of the art is adequate 
in the kinematic regions studied. The strongest caveat to this argument is the possible importance of hadron mass 
effects, which are not so important for tt^, being the lightest hadrons, but which may be relevant for other particles. 
Therefore, for the other particles it may be necessary to account for hadron mass effects in the theory. Furthermore, in 
this connection, it may also be necessary to account for contamination of the sample from decays of unstable particles. 

It is important to note in the discussion above that, due to insufficient information on the systematic effects, even 
a failure to fit certain data points in such a global fit does not necessarily suggest other physics effects, unless the 
theory cannot describe data from different experiments which are consistent with one another. 

Since our previous fits [J, 0] > a- number of measurements have been published by collaborations at RHIC and by 
the CDF collaboration at the Tevatron, which allow for a number of extensions in the knowledge of fragmentation: 
Because the gluon FF only appears at next-to-leading-order (NLO) in the e+e~ + X cross section but at 

leading order (LO) in the pp{p) h"^ + X cross section, inclusion of these data in the purely e+e~ sample would 
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significantly improve the constraints on gluon fragmentation, and may therefore give an FF set suitable for predictions 
of future measurements at e.g. the LHC and RHIC. Furthermore, these charge-sign unidentified measwcements provide 
much needed constraints on the separation between the light quark flavour FFs due to the differences between the 
light quark flavour PDFs. Previously, the only data that could constrain these separations were the OPAL tagging 
probabilities since these are the only data for which light quarks are separately tagged. (Note that quark flavour 
untagged measurements from e+e^ reactions cannot distinguish between quarks with similar electroweak couplings.) 
These data are extracted from single and double inclusive production measurements together with some reliable 
theoretical assumptions, the strongest being SU(2) isospin symmetry between the u and d quark flavours and the 
standard model predictions for the branching fractions of the .Z-boson to each light quark flavour. As required for 
any physical observable, the theoretical definition of a quark flavour tagged measurement is trivially QCD scheme and 
scale independent since it is obtained by setting to zero the electroweak couplings of all quark flavours except that of 
the tagged quark flavour. Finally, valence quark FFs can be constrained by the difference between the production of 
a given species of charged hadron of one sign and the other, the charge-sign asymmetry, from pp reactions at RHIC. 
Since these data depend only on the valence quark FFs, in contrast to the charge-sign unidentified production in e+e^ 
reactions which depend only on the gluon, quark singlet and quark non-singlet FFs, fits to and predictions for the 
charge-sign asymmetry and charge-sign unidentified production can be treated orthogonally. Since the total sample of 
the former data is much lower in quality than that of the latter, the independence of the valence quark FFs from the 
other FFs in this sense and also in the sense that it docs not mix with the other FFs on evolution should be reflected 
in the parameterization used for the FFs at the initial factorization scale. Valence quark FFs will be useful for making 
predictions for other charge-sign separated data such as that at ep colliders — no such data at high negative photon 
virtuality exists at present, however such an extraction from HERA measurements is planned for the future. Such 
data in turn will be useful for improving the currently rather poor constraints on the valence quark FFs. We note 
that there exists data from EMC for charged particle production However, the contamination by charged particles 
other than the ones we are interested in is unknown. Furthermore, there exists data from HERMES for 7r+, 7r~, 
and K~ production from ep reactions but for Q ^ 2 GeV where the validity of the fixed order (FO) approach 
comes into question. However, since such processes are physically very similar to tt^ and production in e~^e~ 
reactions at similar values for ^/s, inclusion of future measurements of the latter type at low ^/s may help to average 
out these unknown theoretical systematic effects. 

In addition, theoretical developments in the calculation of inclusive production observables have occurred since 
the analyses in Refs. [H, Small x divergences can be resumnied within the framework of DGLAP evolution 
which is crucial for improving the description at small x. Since cross section measurements at small x depend on the 
FFs at all larger x values, the inclusion of these measurements will also lead to improved constraints on the FFs at 
the X values currently determined in global fits. Unfortunately, while this procedure to any order is simple, explicit 
results for a full small x resummed NLO calculation do not exist yet. However, effects of the produced hadron's mass, 
which must be treated first, can be incorporated. As we will see later, this allows for e+e" data at smaller x and 
lower centre-of-mass (cm.) energy ^/s to be added to the data sample to be fitted to, and the results convincingly 
demonstrate both the accuracy of the fixed order approach and, particularly for baryons, the fact that hadron mass 
effects are the most relevant small x, low ^/s effects. 

In this paper we repeat the AKK fits of FFs for tt*, K^, p/p Kg and A/A 4] to inclusive production mea- 
surements for these particles. Those fits were intended to be conservative, in that we used only data for which the 
corresponding calculations are reliable. In the fits of this paper we make a number of improvements. Concerning the 
experimental information, we include all e~^e~ reaction data measured below the Z pole. Of these, only TPC data 
were used in Ref. [l| in order to constrain as{Mz). The remaining data were excluded, due to unknown deviations 
from the FO approach, and because the accuracy and number of data points was low. The former problem is handled 
by including hadron mass effects in our calculations and allowing the hadron mass to be fitted. Note that if x and/or 
^/s were too low for the FO approach alone to be valid, this would also have the effect of subtracting out the small x 
and/or low ^/s deviations, such as higher twist. To meet the latter problem we exploit the data further than would 
be done if the normalization error were to be treated as a statistical error added in quadrature or as a normalization 
factor, as is usually the case, by incorporating it instead as a systematic error in a covariance matrix. Note that 
small X, low y/s data impose more constraints on the gluon FF than the larger x, higher y/s ones do. We also resum 
large x logarithms in the perturbative series for the quark coefficient function for e~^e~ reactions, and in the DGLAP 
evolution of the FFs, since this affects all x values in principle and, as we will see later, leads to an improvement in 
the fit over the unresummed case. Furthermore, this modification is simple to apply to our Mellin space calculations. 
The next important update to our fits is the inclusion of charge-sign unidentified hadron production data from pp(j)) 
reactions at RHIC and the Tevatron to further constrain the gluon FF at large x. Similarly as for the e+e~ reaction 
data, we incorporate the systematic normalization errors of all RHIC data via a correlation matrix. The inclusion 
of these data also imposes further tests on universality, namely between pp, pp and, because of the sizable gluon FF 
dependence, e+e" reaction data at smaller x, lower ^/s. 
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In Ref. only the charge-sign unidentified FFs were determined since e+e" reactions are charge conjugate 
invariant, while the valence quark FFs were left completely unknown. Therefore, it is not possible to obtain any 
information whatsoever on the charge-sign asymmetry in fragmentation processes from the AKK FF sets, which are 
independent degrees of freedom. Therefore, the most likely explanation for the discrepancy between the charge-sign 
asymmetry in the BRAHMS data and the corresponding calculation obtained using charge-sign unidentified FFs 
together with certain assumptions relating them to the valence quark FFs is that there is a problem with these 
assumptions. To resolve this issue, we perform a phenomenological extraction of the valence quark FFs from the 
charge-sign asymmetry in the BRAHMS and STAR data, independently of the charge-sign unidentified fits. 

For pp reactions, a charge-sign asymmetry in the produced hadron sample will be observed if the fragmentations 
from the initial protons' valence quarks were to dominate sufficiently over that from the proton's sea partons (gluons 
and sea quarks), since the sea is charge conjugation invariant. Furthermore, the contribution from the fragmentation 
to a particle of a given charge-sign from each of that particle's valence quarks must be positive. We therefore study 
the relative contributions to the charge-sign unidentified charged particle production from the fragmentation of the 
initial protons' valence quarks and sea partons, and to the charge-sign asymmetry from the fragmentation of the 
produced hadron's valence quarks, and look for deviations from our expectations which would signal a problem in the 
fit. 

It is clear that any global fit is incomplete without a full error analysis, including correlation effects, on the fitted 
FFs. Since we wish to develop both the theory and develop and perform the technique for doing this in detail, we 
postpone these analyses to a future publication, while the present paper will be dedicated to the issues given above. 

Since the last AKK analysis, two other analyses have been published in Rcfs. and [ll|. The analyses of this 
paper differ in that we incorporate hadron mass effects to improve the small x and low ^/s regions, resum large 
X logarithms, perform a complete treatment of the normalization errors on the experimental data via a correlation 
matrix, and use pp{p) reaction data including that for Kq and A/A. However, we avoid the HERMES data from 
ep reactions for tt^ and fl] for the reasons given earlier. (Data from pp{p) reactions with transverse momentum 
values as low as ~ 2 GeV are used, however the particle energy is somewhat higher because of the non-zero average 
rapidity, the measurements are for cm. energy > 200 GeV, and the cross section only depends on the FFs at large 
X.) For these reasons our analyses provide a complement to the other analyses. 

The paper is organized as follows. In Section |TT] we explain our choice of data sets used to constrain the FFs. In 
Section Hill we describe the theoretical input for our calculations, and derive new tools for our analyses. Note that 
Sections [Tl] and IIIII are intended to complement rather than repeat the discussions in Refs. 0, 01 ■ AH results are 
presented in Section [TVl and finally the work is summarized in Section IVl 

II. TREATMENT OF EXPERIMENTAL INPUT 

In this section we discuss the motivation for and application of the data used to constrain the FF sets. For further 
details concerning the e+e" reaction data, we refer the reader to Refs. 0, [J- 

A. Experimental data 

The data used for the extraction of FFs for n"^, and p/p particles in Ref. [![ and of FFs for Kg and A/ A 
particles in Ref. ^ were limited to measurements from e"'"e~ reactions for which Xp > 0.1, where Xp is the fraction 
of available spatial momentum taken away by the produced hadron, given in terms of the hadron's momentum 
in the cm. frame and the cm. energy ^/s by Xp = 2\pfi\/\/s. In the analysis of this paper, we make a number of 
extensions to the experimental input. Firstly, to our sample of data from e^e~ reactions we add all available data for 
which ^/s < Mz- Such lower ^/s data suffer from larger theoretical errors and were not considered in Ref. J]] (apart 
from data from the TPC Collaboration, which were necessary for a determination of as{Mz))- For this reason, we 
account for the hadron mass effects 0, which requires only a simple modification to the calculation. Furthermore, 
these older data have larger errors. However, in this analysis we exploit the given statistical information concerning 
the normalization, i.e. the normalization error is treated as a systematic error by incorporating it in a correlation 
matrix, instead of as a random error by adding it in quadrature to the statistical error as is usually done. The general 
treatment of systematic errors and its justification is outlined in Appendix VK\ 

We also include the pp{p) reaction data from inclusive production measurements: in pp reactions at -^/s — 200 GeV 
from the BRAHMS ,9], PHENIX aand STAR collaborations at RHIC, and in pp reactions at = 630 

GeV from the CDF collaboration [IJ] at the Tevatron (the data at ^/s = 1800 GeV are excluded due to their small 
effective x — 2ptI \fs). We impose a lower bound pT > 2 GeV to avoid large theoretical errors at small p^- Any 
errors, systematic and statistical, are added in quadrature, except for the normalization errors which are treated, as 
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for the e+e reaction data discussed above, as systematic errors by using a correlation matrix. For each of Kg and 
A/ A, the CDF data are normahzed by the unknown total cross section, which we therefore leave as a free parameter in 
the fit. Therefore, the CDF data only constrain the shape and relative normalizations of the FFs but not the overall 
normalization. 

We include tt'^ production data from the STAR collaboration [13] in the sample of data for constraining tt^ FFs 
(note that, as in the previous AKK fits, each of our charge-sign unidentified FFs for charged particles and A/A 
is defined to be the FF of a parton of given species and charge for the particle in question added to that for its 
antiparticle), and K'^ production data from STAR [3 and CDF [ll| and K"^ production data from BRAHMS |9|] in 
the sample of data used to constrain both and Kg FFs, by respectively imposing the relations 

Df{x,M^) = ^Dfix,Mf) (1) 

and 

Dfhx,M'f) = ^Df^ix,Mj), (2) 

where j = u, d it i — d, u, otherwise i = j, which follow from the highly reliable assumption of SU(2) isospin symmetry 
between u and d quarks. Due to the nature of the DGLAP evolution, these constraints are independent of Mj, as 
physical constraints should be. 

For each light charged hadron species, we perform two mutually exclusive fits, one fit to the charge-sign unidentified 
data (from e+e" and pp{p) reactions) to constrain the charge-sign unidentified gluon, u, d, s, c and b quark FFs, and 
the other fit to the charge-sign asymmetries from pp reactions at RHIC to constrain the valence quark FFs. The only 
measurements distinguishing A from A that exist come from STAR Ts'] , which are too inaccurate for a reasonable fit 
of valence quark FFs for A/ A, while there are no valence quarks in A'^ and its orthogonal state K^ is not measured. 

For the readers' convenience, we have listed all data sets used in our analyses and their properties in Tables IVljXIl 
to the left of the vertical line. 



III. THEORETICAL FORMALISMS USED FOR THE CALCULATIONS 



In this Section we outline some known theoretical tools that we use for our calculations, not already discussed in 
Ref. [l| , and derive some improvements to these tools. 



A. Production in hadron-hadron reactions 



We briefly outline the theoretical calculation that is used for the pp{p) reactions. A generic reaction hih2 — > h + X, 
where is a hadron, can be described by the quantity 

F,,,.Ax,y,s)^sE^^ipr,y.s) = s ^^^iPT,y,s) (3) 

(exploiting azimuthal symmetry in the second line), where ^/s is the cm. energy, E and p the energy and spatial 
momentum respectively of the detected particle /i, pr its transverse momentum (relative to the spatial momenta of 
hi and h2, which are antiparallel) and y its rapidity, given by 

.-^Inf^, (4) 
2 E -PL 

where is the longitudinal momentum (in the direction of the spatial momentum of hi) of h. The dimensionless 
variable %, which will be convenient for later discussions, is given by 

X=l-V + VW, (5) 

where V and W are the variables typically used in perturbative calculations, related to the usual Mandclstam variables 
S, T and U oi h through 

T 

U 
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Working now in the cm. frame and assuming that the mass of h can be neglected so that E = px coshy, we find that 

X = {2pt/V^) coshy = 2E/^ (7) 

is the fraction of available energy or momentum taken away by the hadron and is therefore the scaling variable of the 
factorization theorem. 

In the factorization theorem at leading twist, the following picture emerges which allows the process to be partially 
calculated perturbatively: In any frame related to the cm. frame via a boost (anti-)parallel to the beam direction 
(for massless hadrons the precise choice of frame is irrelevant in the study below — later when we treat the mass 
effects of the produced hadron we will need to specify a frame), a parton in each initial state hadron hk, k = 1,2, 
moving parallel to it and carrying away a momentum fraction Xk, interacts with the other. The interaction results in 
the inclusive production of a third parton, which subsequently fragments to a hadron moving in the same direction 
and carrying away a fraction x of the parton's momentum. This fragmentation is expressed by writing 

FU.(x,y,'^) = E / ^^-^'U. (f '2/'^'^^/) A' {^,Mj) . (8) 

i ^ 

The equivalent partonic production cross section has been denoted Ff^^j^^- Note that the partonic rapidity is the same 
as the hadronic rapidity, since for massless hadrons y can be approximated by the pseudorapidity 



77 = - In (^tan - j , (9) 

where 6 is the angle which both the produced hadron and the massless fragmenting parton make with the beam in 
the cm. frame. 

Due to the smallness of F^^-^y^^ at sufficiently large xl^^ the FFs at small x do not contribute significantly to Fl^_^f^^. 
To see that (large x) gluon fragmentation contributes significantly to -F^^^^i^ ' ^^^^ measurements of such observables 
are ideal for improving the constraints on gluon fragmentation, consider the dependence of the -F)j^/j, on the parton 
distribution functions (PDFs) /^,, where i labels the parton species and hj — /ii,/i2 labels the initial hadrons, 



Fl,^^ {x,y,s,Mj)=J2 t dx, C dx^F^,^ f^,y,^,^,^^s,Mj) fl\{x^, M])fll{x2, Mj). 

■ ■ J X J — \X\X2 X\ J 



(10) 



The ^7ii2' which are perturbatively calculable [isj, are the equivalent partonic cross sections when the hadrons h, hi 
and /i2 in Fj^^^^ are replaced by the partons i, ii and i2 respectively, with cm. energy squared xiX2S. For simplicity, 
the same factorization scale Mf for all 3 partons has been used. Since the gluon PDF dominates over the quark 

PDFs at small Xi and since F^g at LO is only non zero when i = g, we expect Eq. (fTO|) to be largest when i = g. 

We use the Mellin transform approach as discussed in detail in Ref. [l[ for evaluating the cross sections for the e+e^ 
reactions, since this is fast and numerically accurate. For these reasons we evaluate the pp{p) reaction cross section in 
Mellin space. This can be found by rewriting Eq. ([5]) as a proper convolution (using an obvious simplified notation), 

F{x)= f -f(^)xD{x), (11) 
X \x/ 

so that the Mellin transform, defined by 

f{N)= r dxx'^-^fix) (12) 







for any function f{x), of Eq. (fTTj) is the simple product of the Mellin transform of F{x) and xD{x), 

F{N) = F{N)D{N + 1). (13) 
Note from Eq. pTjl that F{x) is defined only in the range 



X < 2: < 1. 



(14) 
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It will be chosen to vanish outside this range, which is reached for example by the integration region of the Mellin 
transform, Eq. The Mellin transform is inverted according to 



1 



dNx-'^fiN), 



(15) 



c 



where C is any contour which lies to the right of all poles in the complex N plane and which extends to lm(N) — ±oo, 
so that F{x) can be obtained via 



Fix) 



2TTi 



dNx^ F{N)D{N + 1). 



(16) 



Since an analytical determination of the Mellin transform F{N) of F{x) itself cannot be determined (one reason being 
that the evolved PDFs on which it depends may be extracted numerically from, e.g., a grid of values in x and M'j), it 

is necessary to obtain an approximate form for F{x) whose analytic Mellin transform can be obtained. The expansion 
of a function as a weighted sum of Chebyshev polynomials of the first kind [3], whose weights are determined by 
demanding that the expansion and the function agree exactly at the zeros of the Chebyshev polynomials, proves to 
be the most accurate representation of a continuous function as a polynomial of a given degree. In particular, it is 
close to the much harder to calculate minimax polynomial expansion, the polynomial with the smallest maximum 
deviation from the function which is being approximated. In detail, any continuous function f{w) defined over the 
interval — 1 < w < 1 can be approximated by 



/H 



M 



^CkTk-i{w) 



.k=l 



Cl 

2 



(17) 



where M — 1 is the degree of the polynomial, the coefficients c„ are given by 



M 



M 



vr (fc - 



fc=i 



M 



AI 



and the Chebyshev polynomials of the first kind, Tn{w), are defined by 

Toiw) =1 
Ti{w) —w 



(18) 



(19) 



That Eq. (fT7| is exact at the zeroes of the Afth Chebyshev polynomial, i.e. for w values for which Tm{w) — 0, may 
be verified by directly substituting Eq. (fT^ into Eq. pT)) and using the property Tn{cos9) — cos{n9). 

It is not possible to expand F{x) as a Chebyshev polynomial directly, due to the limiting discontinuous behaviour 



lim F{x) cx In ( 1 



(20) 



so we instead apply the expansion to F{x)/ In (1 — x/^) 



F{x) = In (1 



■ M 

E 

.fc=i 



CfcTfc, 



2x-{\ + x) 



(21) 



where we have used the linear map from the Chebyshev polynomial argument ui to x in the range of Eq. (|T4] 

x=2{^'X)w + -{l + x)- 



(22) 



In order to be able to apply the Mellin transform to the right hand side of Eq. (|2ip . the weighted sum over Chebyshev 
polynomials must be rewritten as a weighted sum over the x"" , where a — Q, M — 1. However, this expansion in x is 
numerically very sensitive (45j : Although the magnitude of a Chebyshev polynomial is less than one, the coefficients 
of its expansion in powers of w (and therefore x) are typically much greater than one, growing in order of magnitude 
with M . Therefore, in our calculations, for which M typically needs to be as large as 35 for the data points with x ^ 
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small as ~ 0.02, delicate cancellations among large numbers occur which require somewhat more accuracy than that 
of the 8-byte ("double") precision provided by most computer systems. Various packages for increasing the precision 
of numbers exist; we found the DDFUN90 [4y| package sufficient to achieve both a sufficiently reasonable speed for 
fast fitting and an accuracy to a few parts per mil. The McUin transform of each term In (f — x/^^) 2;" on the right 
hand side of Eq. ([2T|) is calculated as 



u A n— 1 

with the shorthand N' = N + a. The series 

00 „ 

^(^'^) = E^^I A (24) 

n(jy — n) 

is quickly converging for x < I. Si {N) is given by 



5i(iV-f) = ^i=ln(Ar)+7ij--^--^ + -^--^ (25) 



k=l ^ ^ 

where the second equality is the analytic continuation to complex N whose magnitude is large. If |iV| is not large, 
this approximation is still applicable by repeated use of the result 

5i(-7V-l) = 5i(-iV) + l, (26) 

which follows from the first equality in Eq. ([25|) . This last result also allows the Si{—{N + a)) for all the values of a 
to be evaluated quickly from Si{—N). 

Poles created by the denominator in the third line of Eq. (|23p at = n — a for n — I, ...,00 would render the 
inverse Mellin transform inapplicable, however inspection of the numerator reveals that these poles are canceled. 

In the cases where the rapidity y is integrated over, the above procedure can be applied by replacing the value for 
y used in x by the minimum bound on y. 

B. Hadron mass effects 

In this subsection we show how to incorporate the effects of the mass rrih of the produced hadron h into inclusive 
hadron production in e^e^ and hadron-hadron reactions and into their relevant kinematic variables. 

1. e^e~ reactions 

We repeat the study of Ref. [8], including now some discussion on the produced hadrons' energy and on cross sections 
averaged over a range of the scaling variable. The theoretical calculation, which follows from the factorization theorem, 
takes the form 

where z is the light cone momentum fraction and ^ the light cone scaling variable. Specifically, in the cm. frame in 
which we will work from now on, using light cone coordinates V = {V+ = {V° + V^)/V2,V- = {V° -V^)/V2,Vt = 
{V-^, V^)), these variables are defined to be the ratio of the produced hadron's "+" component of light cone momentum 
to the partons' and to the intermediate vector boson's respectively. The energy and momentum fractions measured 
in experiments, given by xe = '^Efi/^/s and Xp = 2\ph\/\/s respectively, are determined from the momentum of the 
intermediate vector boson. 
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of the produced parton. 



0,0 , (29) 

z 



V.= i^-4^^.^]- (30) 



and of the produced hadron, 

Note that the hadron mass effects are accounted for simply by introducing a non-zero "— " component into the hadron's 
momentum. From these resuhs we get immediately our desired relations, 

) ,{ (31) 

A cross section which is measured at a fixed value for some scaling variable x and which is differential in some other 
scaling variable x is calculated from the theoretical result da/d^($^, s) using 

da d^ \d<J , , 

dx dx d^ 

Averaged cross sections are calculated according to 

da\ 1 r^Cs^s) J riCs\s) 

dx/xO<x<x^ x(x , s) ~ x{x" , s) J ^s) dx x{x , s) - x{x" , s) J ^(jio ..-^ 



2. Hadron-hadron reactions 



To incorporate hadron mass effects in pp(p) reactions, we start with the general result from the factorization 
theorem, 

rf</..(P, V^) ^Y. j '^'^^j dx2ft:ixuMj)fl^^^X2,Mj) J dxD^{x,M])dal,^{l,^sl (34) 

where dal_^^^ is the equivalent partonic differential cross section for the production of the massless parton i (which 
will fragment to the produced hadron h), with only the two partons «i and in the initial state. We now work in 
the partonic cm. frame where 1 is parallel to p, since this leads to various simplifications: Firstly, 

. - (35) 
X- ;o + |i| ' 

being the ratio of the produced hadron's + component of light cone momentum to the parton's when the 3-axis is 
aligned with their spatial momenta. Eliminating P = |1| and = ^|pp + gives 



2x|1| = |p| + ^|p|2+to2. (36) 

Secondly, the relation between the experimentally measured observable Ed^(T^_^f^^/dp'^ from the equivalent partonic 
"observable" |l|(i^cr*j,/(iF via Eq. (p4|) can be obtained from the simultaneous results 



E ^|p|2+^^ 



^ (37) 
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where is the sohd angle. Therefore, in the presence of hadron mass, 



3^/1 



E- 



dp3 



J2jdx^J d^2/^(xi,Af|)/^(x2,M|) J dxD^{x,M])\l\ 

ijk 



1 



where the divisor 



i? = 1 - 



(38) 



(39) 



in the integrand is the only modification required to incorporate hadron mass effects. However, we must still obtain 
the relation between |p| of the partonic cm. frame and the quantities px and y of the cm. frame used to define the 
kinematics of the produced hadron in experiments. Since pr has the same value in both frames, the dependence of 
IpI on is easily found to be 



IpI = rriT coshy', 



where 



= -JPt + ml 



(40) 
(41) 



and y' is the rapidity in the partonic cm. frame. Then 

/ 

R = 



1 



V 



The relation between y' and y is 



niT cosh y' + \ rn^ cosh y' — 



(42) 



(43) 



where (jj is the boost rapidity between the partonic cm. and lab frames. To evaluate (jj, note that, in the partonic 
cm. frame, xi,2 are each replaced by a common momentum fraction xie'^ — X2e~'^, which implies 



= In. 



(44) 



Finally, the variables V and W defined in Eq. ^ which are used to define the kinematics of the produced parton 
in the theoretical calculation [43| are obtained from px and y by 



V=l 



W 



(45) 



s 



For completeness, we note from Eq. ([5|) that the above results imply that 

2tot 2m? 

X = ^^coshy 

Vs s 

The relation between y and the pseudorapidity i] is then 

y = sinh^^ ( sinh rj 
For experimental data where a range of 77 values is specified, we use the operator 



1 



V2 - Vi J^i 



drj 



1 



rriT cosh y 

dy- 

\ I m'^ cosh^ y — m'f^ 



(46) 



(47) 



(48) 
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C. Large x resummation 



When the transferred momentum (or energy) fraction x, being equal to Xp for e+e^ reactions and x forpp(p), is large, 
the accuracy of the FO perturbative calculation for the hard part of a cross section and the DGLAP evolution in the 
soft parts is worsened by unresummed divergences occurring in the formal limit a; ^ 1. The accuracy and abundance 
of data may be sufficiently limited at large x that the resummation of these divergences is unnecessary. However, since 
the cross section depends on the hard part in the range x < z < I, where z is the partonic momentum fraction, the 
true "large x" region, being the region where resummation could make a difference, may include significantly lower 
values of x than just those for which x ~ 1. In addition, resummation should reduce theoretical errors. These effects 
of resummation should become more significant at lower y^. Since we use e~^e~ reaction data for ^/s as low as ~ 10 
GeV, and since large x resummation results are easily implemented in FO calculations in Mellin space, we implement 
this resummation in our fits. 

In general, in the series expansion of the hard part W{as,x) of some cross section, these large x divergences take 



the form aj[ln"-' (l- 
take the form a" In"^ 
form 



x)/{l — x)] + , where r = 0, n labels the class of divergence. In Mellin space these divergences 
N. These divergences may be factored out, which results in the calculation of W taking the 



W{a,,N) = W,Uas,N) 



(49) 



where the FO series in parenthesis on the right hand side is free of these divergences since they are all contained in 
Wrcs- W at large N is approximated by Wres when the divergences in >Vres are resummed, which involves writing 
Wrcs as an exponential and expanding the exponent in keeping as In A'^ fixed. We will resum the divergences in 
the quark coefhcient function Cq for inclusive quark production in e~^e~ reactions in the manner of Eq. (j49p . once we 
have obtained to all orders its leading (class r — 0) and next-to-leading (class r = 1) divergences from the results of 
Ref. [13]. The result, Eq. (|55p below, is used in our calculations instead of the FO expression for Cq. 

As we will see, resummation of large x divergences in the e~^e~ reactions typically improves fits of fragmentation 
to charge-sign unidentified hadrons. Formal resummed results also exist for the hard parts of pp{p) reactions inte- 
grated over all rapidity values [48j . However, while the generalization to a given rapidity range can be determined 
approximately, no formal results exist at the time of writing. Therefore, it is not possible at present to perform large 
X resummation in the hard parts of pp{p) reactions. 

The divergences of the e~^e~ reactions (in the MS scheme in which we are working) can be reorganized according 
to the formula [47] 



\nCqiN,asis))= / dz — 

.In 1 



- 1 



^A{as{.s)) + B{as{{l-z).s)) 



0(1), 



(50) 



where 



(51) 



The perturbative series for the quark coefficient function to NLO contains a leading (class r — 0) and next-to-leading 
(class r = 1) divergence, which will be replaced just now by the sum of all divergences belonging to these two classes 
from every order. All of these divergences can be explicitly extracted from Eq. ([50]) by taking 



^(1' =2Cf 
3 



(52) 



B 



(1) 



C 



F- 



Since Eq. ([50)1 is algebraically similar to the resummed quark coefficient function of deeply inelastic scattering [4£ 
we may obtain the divergences of classes r = 0, 1 directly from the MS result in Ref. (soj . 



InC, 



[N. as) 



AW 



A,)ln(l- A,) + A,] + 



^-\n\\-\s)- 



/3o 
A^Pi 
PI 



PI 



PI 



ln(l - A,) 



(53) 
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where = as/3o In A'^. The constant /Sq is that which appears in the expansion of the evolution of as, 

^j[^=/3(a.(M^)) = -f;/3„ar2. (54) 
According to the general form of Eq. ((151), the resulting resummed quark coefficient function that we seek is finally 



i^r— 0,1 

q — 



(l + a.(CW-ar°^iW)). (55) 



The presence in the above equation of Cq "'^ the coefficient of the 0{as) term in the expansion of CJ'""'^ in a^, 
given by 

C-=0-i (1) ^ ^ ^ + (a(i)7b - InTV, (56) 

ensures that the original NLO result is obtained when the whole of the right hand side of Eq. (|55p is expanded in as, 
i.e. when the resummation is "undone" , and therefore also prevents double counting of the divergences. Note that 
there are an infinite number of other schemes which are consistent with this criteria and give the large N behaviour 
of Eq. ((53)) . a typical feature of perturbation theory. 

Equation ([53|) contains a Landau pole when = 1, for which N is real and ^ 1. However, in the inverse Mellin 
transform it is not necessary for the contour in the complex N plane to run to the right of this pole, as it should 
for the other poles, because it is unphysical, created by the ambiguity of the asymptotic series, which in x space is 
essentially a higher twist contribution and therefore negligible. We will keep the contour to the left of this pole, which 
is known as the minimal prescription (5lj |. since this is the most efficient choice for our numerical evaluation of the 
inverse Mellin transform. 

The DGLAP evolution of the FFs also contains large x logarithms, which must also be resummed. We do this 
according to the approach of Ref. [l^l, which uses the fact that the FO splitting functions are already resummed 
[sol [53I I to perform the resummation in the Mellin space evolution analytically. 



D. Evolution with heavy quarks 



NLO matching conditions at the heavy quark flavour thresholds for FF evolution have been derived [54], and occur 
for the gluon FF (if the thresholds differ from the heavy quark pole masses ruq) and the extrinsic heavy quark FF (for all 
choices of threshold). "Extrinsic" here means that the fragmentation from the heavy quark proceeds via perturbative 
fragmentation of the heavy quark to a gluon, followed by non perturbative fragmentation of the gluon to the detected 
hadron. Following the results of spacelike factorization, the intrinsic (non perturbative) component of each heavy 
quark FF, which must be added to the extrinsic one, may be negligible since it is expected to be of ©((Aqcd/w,)^), 
where p is a positive integer. However, the intrinsic component may be important at large x [ssf, and therefore 
we will not neglect intrinsic heavy quark fragmentation here. Since this intrinsic component is unknown, we simply 
parameterize the whole heavy quark FF at threshold and neglect its contribution to the calculation below threshold, 
as in the previous AKK studies. A study of the relative contributions of the extrinsic and intrinsic components of 
heavy quark fragmentation would be interesting, but is not relevant to the goals of this paper, and will therefore be 
left to a future publication. To ensure a continuous gluon FF at each threshold for simplicity, we set the thresholds to 
the heavy quark pole masses nic — 1.65 GeV and mt = 4.85 GeV [5i2.], as opposed to twice these values as practiced 
in previous fits. 



E. Other theoretical choices 



We use NLO results for the calculations of all data. To account for the initial (anti)protons, we use the CTEQ6.5S0 
PDFs |57i]. Consequently, for consistency, Aqcd for 5 active quark flavours is fixed in our fits to the value at which 
these PDFs are determined, 226 MeV. 

Other than the improvements that have been discussed in this section, the theoretical calculations and choices used 
here for the e"*"e~ reaction data are identical to those used in Refs. [l|, Q, and are detailed therein [b^I- There we 
considered only charge-sign unidentified quantities (cross sections and FFs) for which a particle is not distinguished 
from its antiparticle, which are given the label H = 01 h/h and are given by O'' = +0'' (or O''/'' = O'^+O^"). 
These particles were H = 7r±, K^, p/p (H, K° and A/A ;4], requiring 5 independent fits. Since particles are 
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distinguished from their antiparticles in the BRAHMS and STAR measurements for n^, p/p (BRAHMS and STAR) 
and (BRAHMS), we may now also perform fits for the charge-sign asymmetries H — A^h^, for which quantities 
take the form 



QA^h- ^o^ -O" . (57) 

Therefore, for this paper there will be a total of 8 separate fits. Note that these combinations are not the only 
possibility. For tt^, and p/p we could instead consider quantities for which H = and H — . However, in 
that case the corresponding quantities ^■'^ ^ = /2 + (— )0^'='' /2 are statistically contaminated by the highly 
uncertain quantities 0'^='* , so that the more reliable information provided by the quantities O'^ is lost. Secondly, 
the theoretical calculations for O'^ and O'* receive contributions from the same FFs through the mixing occurring 
in the evolution (see Eq. ([50]) below), so that all data must be fitted simultaneously. The combinations H — 
and H = Ach^ that we use avoid all this entanglement. For comparison, we also perform fits without the large x 
resummation, which therefore leads to a total of 16 independent fits. We allow a total of 11 parton species to produce 
the observed hadron through fragmentation, being the gluon, the quarks i = u, d, s, c and b and their antiquarks. 
The FFs are parameterized at a starting scale Mq — \f2 GeV in the form 

Dt (x, Ml) =Nfx'^" (1 - xff (1 + ct (1 - x)^"\ ^^^^ 

D^^'^^ (x, Ml) ^N^^>^^ x'^f'"^ (1 - xff'"^ . 

This is the strongest non perturbative assumption that we will make. (Eqs. ([T]) and ([2]) and Eq. ((6T|l below have more 
physical justification, namely SU(2) isospin invariance between u and d quarks, than the above parameterization.) 
The term proportional to cf was fixed to zero in the previous AKK analysis, and is found to significantly improve 
the fit to the data when freed. The constraints 



D'/^ {x^MD^D'^ {x,Ml) (59) 



and 



Df^'^^ {x, Ml) = -Df^>^^ [x, Ml) (60) 

are exact in QCD and are applied in the fits. For the charge-sign unidentified combinations and the neutral particles, 
the first constraint implies 6 independent FFs, being those for the gluon and the u, d, s, c and b quarks. The second 
constraint implies that charge-sign asymmetry FFs vanish when g is a unfavoured quark (i.e. a sea quark of the 
produced hadron), since the equality also holds without the minus sign, leaving only the FFs for the favoured quarks 
(i.e. the valence quarks of the produced hadron) for charge-sign asymmetry data: for i ~ u and d when H ~ AcTt^ 
and H — Acp/p^ and for i — u and s when H = AcK^. For tt^ we impose the SU(2) u and d isospin symmetry 
relations 

In principle, this relatively precise assumption could also be applied to use FFs to calculate production, so 
that a fit to and data simultaneously could be performed to obtain a single set of FFs for K^/Kg particles. 
However, the differences in the details of these particles' production mechanisms, including large time scale effects, 
may be too significant. 

The hadron mass appearing in the calculation of the pp{p) reaction data must be assigned before fitting, since the 
precalculation of subsection IIII Al depends on it. Each particle's mass will therefore be fixed to its known value. On 
the other hand, since the hadron mass appearing in the calculation of the e^e~ reaction data is not limited in this 
way, it will be freed in the fit. 



F. Contributions of the partonic fragmentations to pp — > (Ac)/i^ + X 

When we come to present our results, we will also verify whether the relative contributions of the fragmentation 
from valence quarks and sea partons to the pp cross sections is consistent with our expectations, to be given below. 
This will also help in determining the sources of the charge-sign asymmetry. Writing out the indices labeling parton 
and hadron production in Eq. gives 

F" = J2 (62) 



t=-nj 
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This can also be taken as the x space result, by defining the product in this equation to be the x space convolution. 
The cross section for H = can be written 

ph± ^ F''^' D'f + F"^^' of + J2 ^^^t' (63) 

where F''^ = F'^ — F^ and refers to sea quarks. (In this subsection, "quark" means a quark or antiquark of the same 
flavour.) The first two terms give respectively the contribution to the production of from the protons' valence 
u and d quarks. These quarks are the source of the charge-sign asymmetry, to be discussed in more detail around 
Eq. Ij65p below. Since there are more valence u than d quarks in the initial protons, the first term is expected to 
dominate over the second for tt* production, since the u and d quark fragmentations are equal, and even more so for 
p/p production, because then u is larger than d quark fragmentation. For K^, d quark fragmentation is unfavoured, 
and therefore the contribution from the protons' valence d quark is expected to be much smaller than from their 
valence u quark. 

The third term or remainder in Eq. (|63p can be regarded as the contribution from the collective sea of the initial 
protons. This term does not contribute to the charge-sign asymmetry because it is charge conjugation invariant. 
Therefore, the larger the third term is relative to the first two, the smaller the charge-sign asymmetry relative to the 
charge-sign unidentified cross section. For hadrons which have non-zero strangeness, or which are superpositions of 
hadrons with non-zero strangeness, which in our set are and A/A in the first case and Kg in the second, the third 
term is expected to dominate: Here, the favoured s quark and u and/or d quark fragmentation from the abundant sea 
occurs, while in the first two terms only the fragmentation from the protons' valence u and/or d quark contributes 
which necessarily involves the production of a heavier s quark. For tt^ and p/p production, for which u and d quark 
fragmentations are favoured, it is not clear whether the first two terms are more important due to the FFs there 
being an order of magnitude larger than the rest, or the third which accounts for fragmentation from the protons' 
abundant partonic sea. In fact, while the decomposition in Eq. provides a simple and direct method to determine 
the underlying partonic physics processes involved, it is "non-physical" , or the physical meaning is arbitrary, since 
each term is factorization scheme and scale dependent. An unambiguous decomposition into factorization scheme and 
scale independent terms is 

Fpp = (Fpp — Fpp) -|- {Fpp — Fpp) -|- {Fpp + Fpp — -Fpp), (64) 

where now the final state superscript "/i^" is omitted for brevity, while the initial state is explicitly indicated by the 
subscript. (Note that Fpp — Fpp. However, we have not made this replacement in order to emphasize that the final 
states differ, by the interchange tt"*" <-> tt~ .) The charge-sign asymmetry originates from the protons' valence u and 
d quarks, represented by the first and second term respectively, where e.g. " means that all PDFs except that 
for the valence u quark are set to zero. In other words, in Mellin space the first (second) term is proportional to the 
square of the u {d) valence quark PDF. Note that these asymmetry generating terms neither depend on the protons' 
sea partons nor receive contributions from interactions between valence u and d quarks. The third term corresponds 
to the charge-sign symmetric contribution, becoming equal to Fpp when the protons' valence quarks vanish. The 
arguments which apply to Eq. (|63p that were given above apply also for Eq. (|64p . It turns out that, qualitatively, 
the relative sizes of the first two terms in Eq. ([M)) are similar to those in Eq. ([55]) . However, for graphical purposes 
we will study the contributions in Eq. ([63]), since the valence quark terms in Eq. ([64]) are typically a few orders of 
magnitude lower than the third term. In other words, the fragmentation from the initial protons' sea partons always 
dominates, even if the charge-sign asymmetry is very significant. 

The charge-sign asymmetry is determined from the FFs according to 

Both terms are factorization scheme and scale independent. Since D^''" is negative and D^"'^^ positive, e.g. an 
excess of 7r+ over tt^ requires a sufficiently large excess of fragmenting u over u quarks relative to d over d. For 
— p/p, all 4 quantities in Eq. I|65p are positive so that a definite excess of p over p is predicted. Likewise, an 
excess of over K~ is predicted, although in this case the second term vanishes since D^"^ = 0. Note that 
there is no dependence on Df"^^ since the s quark can only be found in the protons' sea. In the limit that the light 

quark masses are equal, the resulting SU(3) invariance would imply that this FF is equal to D^"^ . However, this 
symmetry is rather poor for the low energy processes which make up the bulk of these FFs. (SU(3) invariance would 
be useful for the high energy processes, where the light quark masses can be neglected.) Therefore this symmetry 
should be used with care. 
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IV. RESULTS 

The minimized values for the main fits and the fits in which no large x resummation is used are compared 
in Table [D (No comparison is made for the charge-sign asymmetry fits. For these fits only pp reaction data is 
used, so that the difference between the resummed and unresummed fit results only from the differences in the FF 
evolution, which is not as large as that in the quark coefficient function of e+e~ reactions.) For K^^ p/p and A/A, 
the resummation significantly improves the fit. In this sense one can say that the data are sufficiently "large a;" to 
warrant the resummation. 

TABLE L The minimized values in each of the charge-sign unidentified fits. For comparison, the values for the unre- 
summed fit are shown (under "Unres. fit"). 



H 




Main fit 


Unres. fit 


7r± 


518.7 


519.0 




416.6 


439.4 


p/p 


525.2 


538.0 


Kl 


317.2 


318.7 


A/ A 


273.1 


325.7 



From now on we consider only our main fits. Table |TT] shows the results for the hadron masses that are constrained 
by the e+e~ data. Excellent agreement is found for the baryons, which suggests that hadron mass effects comprise 
almost all the deviation from the standard calculation at the smaller x, lower ^/s values of the data considered. While 
the excess in each of the fitted baryon masses is only about 1%, it is somewhat larger for tt^. This is expected because 
the TT^ sample receives large contributions from the decay of the much heavier p(770) to 7r+ -I-tt", while those baryons 
not produced from direct fragmentation will mostly come from decays of their slightly heavier resonances. While the 
fitted masses for and have the correct order of magnitude, there is clearly a large undershoot. A possible 
explanation is the fact that these particles' production mechanisms are so involved that they cannot be even partly 
accounted for by mass effects of the heavier parent particles alone. For example, a kaon (charged or neutral) resonance 
can decay to a pion and kaon simultaneously. Interestingly, the undershoots of the fitted charged and neutral kaon 
masses are very similar, i.e. both are around 154.3 MeV, which is consistent with SU(2) isospin symmetry between u 
and d quarks after all. 

TABLE II: Fitted particle masses used in the calculation of the hadron production from e~ reactions. For comparison, the 
true particle masses are also shown. 



Particle 


Fitted mass (MeV) 


True mass (MeV) 


7r± 


154.6 


139.6 


A-± 


337.0 


493.7 


p/p 


948.8 


938.3 




343.0 


497.6 


A/A 


1127.0 


1115.7 



The minimized reduced (Xdf) values are shown in the last row of Tables FVllXIII In each case, the total x^p is 
in the range 1-2, indicating an overall good quality of fit. However, the Xdf value for each data set can sometimes be 
large due to unknown systematic effects. For this reason some data sets were excluded in Refs. [H, In the fits of this 
paper we have taken the other extreme and have included all data sets in order to enhance the mutual cancellation of 
these unknown systematic effects. In Ref. ill], OPAL quark tagging probabilities were included to improve the quark 
flavour separation of the FFs, particularly in the light quark sector which no other data could constrain at the time 
those fits were performed. In the fits of this paper, the RHIC data also help to separate out the light quark flavours 
due to the different weighting provided by the PDFs. This may be the main cause of the slightly larger Xdf values for 
the individual light quark tagged OPAL probabilities relative to those of the previous AKK fits. Note that the Xdf 
values of the individual light quark tagged OPAL probabilities for K^^ and A/ A are reasonable, as they were for 
the previous AKK fit. This is not surprising since the OPAL data give the predicted strange quark suppression. For 
TT^ and , the fits to the heavy quark tagging probabilities are unsuccessful, which may be due to the large angle 
gluon emission effects that we discussed in Ref. jlj . 
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Turning now to the known systematic effects, the magnitudes of the fitted Xk values, defined in Appendix \X\ to be 
the shift, upwards or downwards, of the measurements resulting from and relative to the Kth source of systematic 
error, are typically less than around 1-2, being the upper limit for a reasonable fit. The most serious exceptions are 
the Xk values for the p/p data from BRAHMS for which 3.25 < y < 3.35 and for the A/A data from CDF and STAR. 
For example, the global tt* production data imply that the central values from TASSO at ^/s =34 GcV overshoot the 
best fit calculation by 3%. All central values from RHIC and Tevatron (CDF) overshoot the best fit calculation, as do 
most e'^e~ reaction data for which a systematic error is given. This may result from contamination of the measured 
sample in each case by other particles, or something else. 

Finally, for completeness and to clarify our choice of parameterizations, the central values of all parameters are 
listed in Tables IIIII and IIVI Due to the low factorization scale at which these FFs are parameterized, extractions of 
our FFs at arbitrary momentum fraction and factorization scale should be performed using the grids and FORTRAN 
routines referred to in Section |Vl They should not be extracted by using these parameters directly in a NLO evolution 
routine, which will in general differ from our routine by a potentially large NNLO error at Mf — 0{1) GeV. For the 
same reason, any physical interpretation of these parameters should be avoided. However, we note for the charge-sign 
unidentifed FFs in Table IHII that whenever a large value (> 1000) for Ni is obtained, it is compensated for by a 
large value (> 9) for 6^, which signifies a large statistical correlation between these parameters in these cases, and 
consequently poor constraints at large x. Conversely, for the charge-sign asymmetries in Table HVl poor constraints at 
small X were signified in previous fits by a large correlation between the N.^"^ and a^"^ , which is why the a^"^ 
were fixed to zero for these FFs in our final main fits. Note that SU(2) isospin symmetry between u and d was used 

to constrain D^"^ — D^"^^ , while due to poor experimental constraints we fixed D'^"^^^ — O.BDu"^^^ ■ 



TABLE III: Values, to 2 decimal places, of the parameters of the charge-sign unidentified FFs at Mj = \/2 GeV for which the 
values are minimal. We fix (a;,M|) = DZ {x,Mf) as dictated by SU(2) isospin symmetry between u and d quarks, 
and also fix c"* = d^ =0 since they were rather independent of ■ 



Parameter 


^± 




p/p 




A/ A 




247.80 


16.11 


16155.68 


1.64 


26.92 


ag 


1.93 


2.13 


7.26 


0.84 


4.49 


b. 


6.14 


3.28 


9.07 


4.11 


5.18 


Cg 


0.96 


0.78 


2.04 


1.18 


3.58 


dg 


-0.53 


2.26 


-0.43 


-0.07 


-1.31 


Nu 


0.32 


1.66 


0.49 


3781.89 


0.60 


flu 


-2.07 


0.22 


-0.05 


4.68 


-0.27 


bu 


0.96 


3.55 


1.84 


16.79 


2.25 


Cu 


-0.81 


0.50 


-0.24 


2.34 


0.01 


du 


2.91 


-1.74 


-0.01 


-0.26 


-2.67 


Na 


^Nu 


3.10 


0.03 


121.78 


0.71 


dd 


= au 


-0.29 


-2.61 


3.89 


-0.62 


bd 


= bu 


6.71 


0.69 


9.68 


3.32 


Cd 


= Cu 


-0.07 


-0.91 


1.67 


2.56 


dd 


= du 


5.52 


0.46 


0.26 


19.80 


Ns 


152607.12 


0.82 


3574.00 


659.46 


3.65 


as 


7.34 


-0.04 


10.57 


6.31 


0.00 


bs 


12.29 


1.62 


16.87 


6.80 


4.69 


Cs 


(fixed) 


1.16 


39.06 


1.53 


0.18 


ds 


(fixed) 


0.06 


-6.55 


0.45 


-3.67 


iVc 


0.33 


12.06 


43.30 


6.82 


6.68 




-2.05 


0.99 


2.35 


2.19 


0.43 


be 


2.61 


4.77 


9.36 


5.87 


5.29 


Cc 


-0.88 


5.45 


15.04 


0.92 


0.78 


dc 


2.13 


6.52 


13.74 


-0.35 


-0.07 


Nb 


1.25 


15.72 


6.81 


17.23 


35.20 


ttb 


-0.45 


0.96 


0.48 


1.32 


0.60 


bb 


4.37 


7.94 


11.89 


12.17 


18.91 


Cb 


17.48 


21.05 


0.43 


0.84 


1.45 


db 


10.79 


11.38 


0.00 


-0.02 


-5.59 
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TABLE IV: As in Table|IIIl but for the charge-sign asymmetry FFs at Mf = y/2 GeV. We fix D'^"'' (x, Mj) = D^"'' {x, Mf) 
as dictated by SU(2) isospin symmetry between u and d quarks, and also fix af"^ — af^"^^^ = since they were rather 
independent of x^- In addition, we fix D'^"^ = as dictated by charge-sign symmetry. 



Parameter 


A,7r± 




AcP/p 


Nu 


153.26 


0.73 


52301.42 


au 


13.00 


(fixed) 


11.06 


bu 


2.51 


1.45 


6.74 


Nd 


^N.u 


(fixed) 


= 0.5iVu 


ad 


= au 


(fixed) 


= au 


bd 


— bu 


(fixed) 


= bu 



To illustrate the various features of the theoretical ealculations, we inelude a number of figures. First we study pp(j)) 
reactions: For charge-sign unidentified hadrons, we show the comparisons of the calculation with the experimental 
data in Fig.[Tl the effects of hadron mass and the relative theoretical errors in Figs. [Hand [31 and the calculation using 
other FF sets in Figs. [4] and [5] For the charge-sign asymmetries, these features are respectively shown in Fig. [6l in 
Figs. [7| and [5| and in Figs. and [TUl Finally, the contributions to the production from the fragmentations of the sea 
and valence partons of the initial protons at RHIC for charge-sign unidentified hadrons are shown in Figs. [TT] and [T2l 
and for the charge-sign asymmetries in Figs. [T51 and [Til Next we study e+e" reactions: We show the comparisons of 
the calculation with some of the experimental data for a representative range of -^/i in Fig. [151 the effects of hadron 
mass, large x resummation and the relative theoretical errors in Figs. [TBI and [TTl and the calculation using other FF 
sets compared with data at y/s — 91.2 GeV, at y/s = 14 GeV and at various values of -^s in Fig. [THl in Fig. [T^ and in 
Fig. [201 respectively. Finally, the calculations using the various FF sets compared with the OPAL tagging probabilities 
are shown in Figs. [2TH251 . Lastly, we show the various FFs at y/s = 91.2 GeV directly: the charge-sign unidentified 
FFs in Figs. [^511501 and the charge-sign asymmetry FFs in Figs. [5TH551 

For inclusive tt* production in pp{p) collisions (Fig. [T]), the calculation agrees well with the PHENIX and STAR 
data. Although the calculation undershoots the BRAHMS data for both rapidity ranges, this is not so serious when 
the systematic effects are taken into account since the magnitudes of the fitted Xk (see Table [V| are less than 2, as 
discussed earlier. This undershoot is more serious for the p/p fit. The FF sets for each of K'^ and fit the data for 
both i^g and well, which is consistent with SU(2) isospin invariance. The calculation with the A/A FF set agrees 
reasonably well with the CDF data at the larger pT values only, and significantly undershoots the STAR data. As 
expected, hadron mass effects are negligible in tt^ production data regardless of the rapidity (Figs. [Hand [31). A weak 
effect is seen in the slightly heavier /Kg production data, which does not depend on rapidity but which is largest 
at the smaller px values as expected. The suppression due to hadron mass effects is strongest for the heavier baryons 
p/p and A/A, and is slightly larger at higher rapidity. We also perform the calculations using the AKK 0,13, DSS 
|lll |. DSV [5g (for A/A) and HKNS [lC| FF sets (Figs. [H and [5l). In these cases the hadron mass is set to zero, as was 
done in those fits. It is expected that all FFs including the gluon FF for tt^ production are reasonably well constrained 
since the e'^e~ — > tt^ -|- X data are fairly precise and lead to good agreement with pp reaction data. The differences 
in the calculations using different FF sets are then most likely due to theoretical errors from different choices for the 
various (NNLO) ambiguities in the calculations, such as the calculation of as(/i^), which explains why the calculations 
for TT^ at low and high rapidity are similar relative to the theoretical error, and also why the calculation using the 
previous AKK set is close to the one using this paper's FF set. However, for the HKNS calculation deviates 
from the rest at high rapidity. The calculation with the previous AKK set for p/p at high rapidity deviates from the 
others but gives good agreement with the BRAHMS data. However, the other FF sets agree better with the STAR 
data at lower rapidity. For Kg, the description of the data is better with the fit of this paper than with the previous 
AKK one. However, for A/ A, the situation is the opposite, and in fact there is a strong discrepancy between theory 
and data when using this paper's and the DSV FF sets, implying a possible inconsistency between the pp and e'^e~ 
reaction data for A/A production. On the other hand, rather good agreement is obtained when using the previous 
AKK FF set, because the gluon FF for A/A at the initial scale in Ref. was fixed to 1/3 that for the AKK proton 
for this purpose. Ultimately, a determination of the error on this prediction from the experimental errors, including 
correlation effects, on the FFs would better determine whether an inconsistency really exists. 

The BRAHMS data for which 3.25 < y < 3.35 (Fig. [5]) provide most of the constraint on the charge-sign asymmetry 
FFs, while the constraints from the STAR data are rather poor. The description of the A^p/p data from BRAHMS 
for which 3.25 < y < 3.35 is particularly poor. The AcTt^ data are much less precise than the tt^ data, which is 
due to the similar yields of 7r+ and 7r~ relative to the experimental error. The AcK^ and Acp/p data do not suffer 
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this problem as much, particularly the ^cVlv data, since the yields for each charge-sign, particularly in the case of 
p/p, are significantly different. For the higher rapidity data at BRAHMS, the theoretical error is slightly lower for 
AcTT^ (Fig. [7|) than for tt^ (Fig. [2|), and lower to a greater degree for Ac-ftT^ than for , and for IS.cvIv than for 
p/p, which may result from some cancellation of the theoretical error between the cross sections for each charge-sign, 
whose calculations are similar. However, the theoretical errors for AcTt^ (Fig. [8]) and tt^ (Fig. [3]) data at STAR are 
similar, while the theoretical errors for Acp/p are larger than for p/p, which may be due partly to (hidden) theoretical 
errors at low rapidity. Note that, perhaps for the same reason, some cancellation of hadron mass dependence also 
occurs for Acp/p relative to p/p at high and low rapidity. For the higher rapidity data (Fig. ^ the calculations 
with this paper's FF sets are somewhat higher than the others (at higher pr for AcTt^), presumably due to larger 
experimental errors in the data fitted to. This occurs to a much greater degree for A^p/p at low rapidity, probably 
because, in contrast to the other FF set fits, we impose additional constraints on the A^p/p FFs with the BRAHMS 
data for which 2.9 < y < 3, with which the calculations using the other FF sets disagree strongly. The calculations 
for AcTT^ for STAR kinematics using this paper's FF set and the DSS FF sets are similar, and somewhat different to 
the calculations using the HKNS FF set (Fig. [T0|) . Note that in this latter set the charge-sign asymmetry FFs were 
constrained using theoretical assumptions only. 

For TT^ , the fact that charge-sign asymmetry is observed at BRAHMS and not at STAR is explained by the greater 
importance of the charge-sign asymmetric fragmentation from the initial protons' valence u and d quarks (Figs. [TT] and 
[T2l) . The charge-sign asymmetry at STAR should become visible at sufficiently large pr- The excess of fragmentation 
from valence u quarks over valence d implies an excess of tt^ over 7r~, because then the first term in Eq. (|65p . which 
is positive, is larger in magnitude than the second, which is negative (Figs. [T51 and [H]). (If the magnitude of the 
"d -I- d" contribution to the excess of 7r+ over 7r~ were larger than that of the "m -I- -0" , this excess would become 
negative). The lower fractional contribution from the initial protons' valence quark fragmentations to the cross section 
for the STAR data implies that this charge-sign asymmetry increases with rapidity. As expected from the discussion 
following Eq. (j63|) . the contribution to production from fragmentation of valence d quarks is negligible. While 
the production of exceeds that of K~ ^ at low rapidity the small contribution from valence quark fragmentation 
suggests that any charge-sign asymmetry is difficult to observe. For p/p at high and low rapidity, valence u quark 
fragmentation dominates over d, more so than for tt^ as expected. It is clear that the excess of p over p (relative to 
the p/p production) increases with and rapidity. At high rapidity, the valence u and d quarks contribute to the 
excess of p over p as expected. However, at low rapidity the contribution from the fragmentation of valence d quarks 
to p exceeds that to p. This is inconsistent with our expectations. Note that there may be larger hidden theoretical 
errors at lower rapidities. This was also the reason given earlier in this Section for the larger theoretical error for 
Acp/p than for p/p. In addition, the propagated experimental errors on the FFs may still allow for this contribution 
to the charge-sign asymmetry to have the opposite sign. This is suggested by the fact that the cross section is much 
smaller at these lower rapidities. 

We now turn to inclusive production in e+e^ collisions. Good agreement is found with the OPAL data at -^i — 91.2 
GeV for all particles (Fig. [151), ^ indicated by the Xpp values in Tables IVl I Vlli IXII and IXIIl except for p/p, which 
is indicated by the high x^p value in Table Hxl Reasonable agreement with the TASSO data at -^/s = 34 GeV is found 
for all particles. Note for example that, for these data, significant disagreement is found for tt^ at the two largest x 
values, despite the value of \k being positive (see Table |V]) , because in fact the theory would otherwise overshoot 
the rest of these data, which furthermore are more accurate. For all particles the calculation tends to undershoot the 
less accurate TASSO data at \fs = 14 GeV. However, looking at tt* for example, the value for \k for these data is 
strongly negative and its magnitude is not unacceptably large. A fitted mass for p/p results in a strong suppression of 
the cross section for a: ^ 0.5 (Fig. [T6)) . implying that hadron mass effects cannot be neglected in this region. On the 
other hand, hadron mass effects for tt^ are negligible. In general, the effect of large x resummation is to enhance the 
cross section at large x. However, note that there is some suppression at small x. Formally, this suppression is less 
important than the effect of resummation of small x divergences arising from soft gluons. The theoretical error both 
with and without resummation is similar (on a logarithmic scale), although it must be remembered that the same 
choice in both cases for the variation in k might not be appropriate. At large ^/s (Fig. [17]), the mass effects for p/p are 
neghgible as expected. Resummation effects are also reduced at higher ^/s. At = 91.2 GeV, all FF sets give similar 
results at a; < 0.4 (Fig. [T5)) . but there are of course some differences at large x where the experimental uncertainties 
are larger. This paper's fit and the calculation with the HKNS FF set for tt^ prefer the SLD data around x ~ 0.6, 
but, as for the other FF sets, the FF set of this paper tries to fit the DELPHI, OPAL and SLD data at a; ~ 0.8. The 
result for at large x tries to fit both the OPAL and SLD data points at a; ~ 0.8. For p/p, the OPAL data at all 
X values are clearly inconsistent with the others, but this does not affect the good agreement with the other data. 
However, the calculation with this paper's FF set overshoots the others at values of x above those for the data. For 
A/ A, the calculations are all significantly different at large x. Note that the DELPHI data are inconsistent with the 
others at small x. The agreement between the calculations for tt^ at lower ^ = 14 GeV (Fig. [T9|) is similar to that 
at ^/s = 91.2 GeV. Note that data on inclusive tt* production is generally much more accurate and abundant than 
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for other particles. At intermediate x, slightly more disagreement among the calculations exists for for ^/s = 14 
GeV than for ^/s = 91.2 GeV. Strong disagreement among the calculations at smaller x and y/s = 14 GeV for p/p 
is found. However, the calculation with this paper's FF set describes the smaller x HRS and TPC data at y/s = 29 
GeV better (Fig. I^O)) . and also the TOPAZ data at y/s — 58 GeV where the calculation also "points to" the data for 
which X < 0.05. For Kg, both this paper's FF set and the previous AKK FF set give similar results (Fig. [T9)) . The 
calculation with this paper's FF set for A/A are significantly different to the others at the smaller x values, although 
the agreement here with the TASSO data at ^/s — 14 GeV is better. 

Since the OPAL tagging probabilities are physical, they give a better indication than the FFs themselves do of 
how well the fragmentation from the individual quarks is currently understood. (Nevertheless, for completeness we 
will also study the FFs themselves next.) For (Fig. [H]), u and d quark fragmentations are well constrained, and 
are consistent with the corresponding OPAL tagging probabilities, while s and b quark fragmentations are badly 
constrained at large x. The e~^e~ reaction data constrain the fragmentation of these quarks (and the d quark) 
less than the others because of their smaller coupling to the electroweak boson. Since the corresponding tagging 
probabilities have large errors and untagged or light quark tagged data cannot constrain the difference between s and 
d quark fragmentations, the only other constraints, in the case of this paper's and the DSS FF sets, come from the pp 
reaction data, which must therefore be rather poor for the s quark fragmentation where the two calculations disagree 
markedly. However, additional constraints for c and b quark fragmentation are provided by light, c and b quark tagged 
cross section measurements from DELPHI, SLD and TPC (see Table IVT) . The calculations overshoot the measured 
tagging probabilities at the lower x values, particularly the heavy quark ones, which may be due to inconsistencies 
between the theoretical and experimental definitions in this region as discussed in Ref. ;li]. For (Fig.IH]), d quark 
fragmentation is poorly constrained, and unfortunately little more constraint is provided by the tagging probabilities. 
However, the other quark fragmentations are better constrained since their calculations using the different FF sets 
yield similar results. Note that for both tt^ and K^, the favoured u quark fragmentation is rather well constrained, 
while the unfavoured light quark fragmentation (s and d respectively) is not. For p/p (Fig. I23p . more discrepancies 
are found at large x, particularly for the heavy quarks although here the calculations with this paper's and the DSS 
FF sets are relatively less separated. Fragmentation to p/p from u quarks is fairly well understood, and also from the 
other quarks except at the large x values. For Kg (Fig. lMj) . the calculations using this paper's and the previous AKK 
FF sets are similar. Good agreement is obtained with the (most accurate) s and b measured tagging probabilities, 
even at the smallest x. For A/A fFig. [^5]). as for Kg for all quarks, the calculation with this paper's FF set is close 
to that using the previous AKK FF set for s and c quark tagging, but otherwise these and the calculation with the 
DSV FF set are quite different. 

Finally, we compare the FFs at Mf = 91.2 GeV fFigs. [^5H5D)) . In general, the relative behaviour of the FFs at lower 
Mf values is qualitatively similar. For each parton, this paper's FFs are typically lower than the others because the 
large x resummation enhances the cross section. Fragmentations from the gluon are poorly constrained, even those 
from fits where pp reaction data is used (comparing this paper's gluon FFs with the corresponding DSS ones where 
applicable), because although these data depend more strongly on the gluon FF than the e+e" reaction data, the 
former data are much more limited than the latter. For tt^, K^ and Kg, the relative behaviour of the quark FFs 
for each flavour are similar to their tagging probability counterparts discussed above. However, excluding HKNS, for 
p/p there is better agreement at large x among the u and d quark FFs than among their tagging probabilities, the 
uncertainty in the latter case presumably coming from the contributions of the other quark FFs. For A/ A, the FFs 
are rather similar at the lower x values, but, apart from the agreement between this paper's and the previous AKK 
FFs for the s and c quark at large x, are otherwise markedly different. 

For AcTT^ , the similarity between the DSS and HKNS results for the u and d (Fig. [ST]) quark FFs and the u quark 
FF for AcK^ (Fig. EJ) and for A^p/p (Fig. ^ is to be expected since similar non perturbative assumptions were 
made for these fits while for the fits of this paper we have let the data decide. For AcTt^, this paper's quark FFs are 
typically much smaller than the other FFs, while the situation is the opposite for K^ and even more so for A^p/p. 



V. CONCLUSIONS 

We have determined charge-sign unidentified FFs for tt^, K^, p/p, Kg and A/A, and FFs for the charge-sign 
asymmetries A^tt^, A^K^ and A^p/p from e+e^ and pp{p) reaction data. Relative to the previous AKK fits [l|, 0], 
we added e'^e~ reaction data at smaller x and lower ^/s, which provides stronger constraints on the gluon fragmentation 
that enters at NLO. To account for deviations from the FO calculation in this region, we included hadron mass effects 
in the calculation, with the mass in the calculation for the e+e~ reaction data left as a free parameter in the fit. For the 
baryons p/p and A/A we obtained excellent agreement with the true masses. Therefore these particles are probably 
produced almost exclusively from direct partonic fragmentation, and therefore may be key to our understanding of 
the partonic fragmentation process. This is then a strong reason for the inconsistencies between the descriptions of 



19 



e+e" and pp reactions in baryon production to be resolved. In particular, the description of the STAR data for A/A 
is poor, while the contribution from the initial protons' valence d quarks to the charge-sign asymmetry for p/p from 
STAR is negative. However, good agreement with the BRAHMS data for p/p was found in general. For the mesons 
TT^, and Kg we obtain the correct order of magnitude. The overshoot (undershoot) in the case of tt^ (if^ and 
K'g) suggests that a significant portion of the sample originates from the decay of heavier particles (complicated decay 
channels). We implemented large x resummation in e+e^ reactions and in the DGLAP evolution of the FFs, which 
for most particles leads to a significant improvement in the fit. We included RHIC pp reaction data for all particles 
and Tevatron pp reaction data from the CDF collaboration for K'g and A/ A to improve the constraints on the gluon 
fragmentation, the quark flavour separation and also, in the case of the RHIC data, to determine the charge-sign 
asymmetry FFs. Hadron mass effects were included in the calculation of these data as well. For both e+e" and pp{p) 
reaction data, the normalization errors were accounted for by including them in the correlation matrix used for the 
calculation of x^. The corresponding \k values were determined, and were typically in the reasonable range \\k \ ^ 2. 
In order to prevent the large errors of the statistically lower quality charge-sign asymmetry data propagating to all 
FFs, we fitted the valence quark FFs on which they depend separately from the other FF degrees of freedom. This 
avoids the need to make unreliable non perturbative assumptions which will ultimately affect all fitted FFs in unclear 
ways. 

In order to make predictions, our FFs sets over the fitted range 0.05 < z < 1 and Mg < Al'j < 100000 GeV^ can be 
obtained from the FORTRAN routines at http://www.desy.de/~simon/AKK.htinl, which are calculated using cubic 
spline interpolation on a linear grid in z and linear interpolation on a linear grid in InM^. We note a number of 
crucial points relevant to the use of our FF sets. Firstly, it is neither incorrect nor even inconsistent to use NLO FF 
sets such as ours in a LO calculation. The result will simply be to LO only. Secondly, while our mass corrections 
are desirable, it is not necessary to include them in calculations using our FF sets, which would then be at least as 
accurate as cross sections calculated using other FF sets that were extracted without the use of mass corrections. 
Finally, the FF sets here were fitted at a low factorization scale Mq = GeV, where the N^LO error on our NLO 
evolution may be quite sizable. It is therefore incorrect to evolve our FFs using a different NLO evolution procedure 
(e.g. one in which the DGLAP equation is solved numerically). Rather, the FFs supplied at the web site just given 
should be used. 
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APPENDIX A: CORRELATED ERRORS 

In this appendix we assume symmetric errors. In the absence of systematic effects, the probability for observables 
fi to take values between /* and /* + df*, given measured values /? ± CTj, is proportional to exp[— x^/2]d/*, where 

^^=i:(^)'- (Ai) 

The Kth source of systematic uncertainty will cause ff to be shifted to f^ + Xxcrf^ , where the probability density 
in Xk is proportional to exp[— A^/2], which defines each systematic uncertainty af. Therefore, systematic effects 
modify to 

= E ( ^''■^' '^""""^ )' + E ^K- (A2) 
The most likely values of the Xk occur at 



(A3) 
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Solving these equations for the \k gives 
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where the covariance matrix 
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and Eq. (jA2p becomes 
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TABLE V: Summary of the fit for vr* production. Details of the data used are shown to the left of the vertical double fine, fit 
results to the right. The columns labeled "Norm." and "Shift" respectively give the experimental normalization and "fitted" 
(i.e. using the fitted Xk) normafization error(s) and as a percentage. In the column labeled "Properties", "i tagged" means the 
light quarks are tagged. The BRAHMS data have 5 sources of normalization error, whose values are shown for pr above (below) 
3 GeV. The normalization error of 11.7% on the STAR data follows from the total cross section measurement a = 30 ±3. Smb. 



Collaboration 


Properties 


(GeV) 


# 
data 


Norm. 

(%) 


2 

Xdf 




Shift 
(%) 


TASSO fl5] 


untagged 


12 


5 


20 


0.56 


0.26 


5.3 


TASSO fl6] 


untagged 


14 


10 


8.5 


0.91 


-1.24 


-10.5 


TASSO fl6] 


untagged 


22 


1 


6.3 


0.00 


0.05 


0.3 


HRS [17] 


untagged 


29 


6 




1.37 






TPC fl8] 


1 tagged 


29 


9 




0.30 






TPC fl8] 


c tagged 


29 


9 




0.68 






TPC fl8] 


b tagged 


29 


9 




1.23 






TPC fl9] 


untagged 


29 


27 




1.09 






TASSO fl5] 


untagged 

Oo 


30 


4 


20 


0.61 


0.72 


14.5 


TASSO [20] 


untagged 


34 


10 


6 


1.11 


0.61 


3.7 


TASSO [20] 


untagged 

OO 


44 


7 


6 


2.01 


0.67 


4.0 


TOPAZ [21] 


untagged 

Oo 


58 


8 




0.85 






ALEPH [22] 


untagged 

OO 


91.2 


22 


3 


0.57 


-0.55 


-1.6 


DELPHI [23] 


1 taeeed 


91.2 


17 




1.76 






DELPHI [23] 


b tagged 


91.2 


17 




1.76 






DELPHI [23] 


untagged 

Oo 


91.2 


17 




1.76 






OPAL [5] 


u tagged 

oo 


91.2 


5 




5.78 






OPAL [5] 


d tagged 


91.2 


5 




5.43 






OPAL [5] 


s tagged 


91.2 


5 




4.43 






OPAL [5] 


c tagged 


91.2 


5 




7.25 






OPAL [5] 


b tagged 


91.2 


5 




10.92 






OPAL [24] 


untagged 


91.2 


20 




1.26 






SLD [25] 


1 tagged 


91.2 


28 




0.79 






SLD [25] 


c tagged 


91.2 


28 




0.78 






SLD [25] 


b tagged 


91.2 


28 




0.71 






SLD [25] 


untagged 


91.2 


28 




0.57 






DELPHI [26] 


untagged 


189 


3 




4.67 






RR AHMS f9] 


[2.9,3] 


200 


8 


11,7,8(13), 
2,1(3) 


0.90 


-1.67, -1.06, -1.20, -0.30, -0.14 


-18.4, -7.4, -9.6, -0.6, -0.1 


y e [3.25, 3.35] 


7 


2.83 


-1.91, -1.21, -1.81, -0.35, -0.34 


-21.0, -8.5, -14.5, -0.7, -0.3 


PHENIX [2] (tt") 


\v\ < 0.35 


200 


13 


9.7 


0.48 


-0.01 


-0.1 


STAR [12] (tt") 


■n = 3.3 


200 


4 


16 


0.72 


-0.58 


-9.4 


STAR [12] (tt") 


ri = 3.8 


200 


2 


16 


0.54 


-0.27 


-4.4 


STAR [3] 


\y\ <0.5 


200 


10 


11.7 


0.48 


0.08 


0.9 


Total 






382 




1.36 







TABLE VI: As in Table FVl but for the charge-sign asymmetry AcTt*. 



Collaboration 


Properties 


(GeV) 


# 
data 


Norm. 

(%) 


2 

Xdf 


Aa- 


Shift 

(%) 


RR AHMS [9] 


[2.9,3] 


200 


8 


11,7,8(13), 
2,1(3) 


0.66 


-0.59, -0.38, -0.48, -0.11, -0.07 


-6.5, -2.6, -3.8, -0.2, -0.1 


y G [3.25, 3.35J 


7 


0.47 


-0.39, -0.25, -0.25, -0.07, -0.02 


-4.3, -1.7, -2.0, -0.1, 0.0 


STAR [3] 


\y\ < 0.5 


200 


10 


11.7 


0.07 


-0.05 


-0.5 


Total 






25 




0.37 
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TABLE VII: As in Table |Vl but for K"^ production. 



Collaboration 


Properties 


V * 
(GeV) 


41 

IT 

data 


Norm 

(%) 


Xdf 


\k 


Shift 

(%) 


TA QQPi ri 

l/looL^ IIOJ 


1 J 1 

untagged 


1 


Q 

o 


on 

zU 


n 00 
u.yz 


-U.Oo 


-i i.O 


lAOOKJ |1DJ 


untagged 


1 A 




y 


o.O 


0/11 
z.4i 


n 01 

U.zi 


1 s 

i.O 


l/Voo'-J |1D] 


untagged 


00 


7 


D.o 


1 01 
i.Zi 


-i.OO 


Q s 
-y.o 


WPQ n 71 


untagged 


OQ 


7 




1 1 

i.iO 






A/TAT51x"TT [071 
iVi/Vxlxvii 


untagged 


oo 


o 
z 


1 
iz 


C\AA 
U.44 






TPn ri Ql 


untagged 


00 

zy 


ZD 




u.uu 






TA QQr> ri 


untagged 


QO 



z 


on 

zU 


n /IS 

U.4o 


n /IO 

-U.4Z 


-o.O 


T A QQM Fonl 


untagged 


o4 





D 


n on 
U.zU 


-U.oD 


O 

-z.z 


TOP A V [01 1 
IVJJr^AZj IziJ 


untagged 


Oo 







n 1 Q 
U.io 






A T TT'PH Too OSl 


untagged 


01 

yi.z 


1 s 
io 


o 
o 


u.oo 


n "^7 


1 1 
i. i 


PlPT PWT roQi 
UsliLirrll. I^oj 


1 tagged 


01 

yi.z 


1 7 
i / 




n 77 








b tagged 


01 

yi.z 


1 7 
i i 




n 77 






■ntTT PMT FOQl 
Uil/ijJr^rli |zoJ 


untagged 


oi o 
yi.z 


1 7 
i / 




n 77 






OP AT [c;! 


u tagged 


01 

yi.z 







111 
i.ii 






OPAT \K] 
\Jri\Li \o\ 


d tagged 


01 

yi.z 







n 7S 

U. ( o 






OPAT \K] 


s tagged 


01 

yi.z 







1 07 

i.y ( 






OP A T IK] 


c tagged 


oi o 
yi.z 







71 

0. ( i 






OPAL f5] 


b tagged 


91.2 







fi 07 

o.y ( 






OP At To /1 1 
U-tAL |z4J 


untagged 


m o 

yi.z 






U.4 ^ 






SLD [251 




91.2 


28 




1.82 






SLD [25] 


c tagged 


9L2 


28 




0.96 






SLD [25] 


b tagged 


91.2 


28 




1.93 






SLD [25] 


untagged 


91.2 


28 




0.46 






DELPHI [26] 


untagged 


189 


3 




1.10 






RRAHMS m 


y G [2.9, 3] 


200 


8 


11,7,8(13), 
2,1(3) 


1.14 


-1.25, -0.80, -0.35, -0.23, 0.00 


-13.8, -5.6, -2.8, -0.5, 0.0 


y G [3.25,3.35J 


6 


2.93 


-0.78, -0.49, -1.45, -0.14, -0.25 


-8.5, -3.5, -11.6, -0.3, -0.2 


CDF [14] (Ag^ 


\v\ < 1 


630 


37 


10 


0.59 


-2.14 


-21.4 


STAR [3] {K'^) 


\y\ < 0.5 


200 


9 


11.7 


1.04 


-1.70 


-19.8 


Total 






346 




1.20 







TABLE VIII: As in Table IVlIl but for the charge-sign asymmetry Ac/f^ 



Collaboration 


Properties 


(GeV) 


# 
data 


Norm. 

(%) 


Xdf 


Aa- 


Shift 

(%) 


BRAHMS f9] 


y G [2.9, 31 


200 


8 


11,7,8(13), 
2,1(3) 


1.15 


-0.59, -0.38, 0.03, -0.11, 0.04 


-6.5, -2.6, 0.2, -0.2, 0.0 


y G [3.25, 3.35] 


6 


1.07 


-1.13, -0.72, -1.02, -0.21, -0.14 


-12.4, -5.0, -8.1, -0.4, -0.1 


Total 






14 




1.12 
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TABLE IX: As in Table |Vl but for the fit for p/p production. 



Collaboration 


Properties 


V 

(GeV) 


ir 
data 


Norm. 

(%) ' 


Xdf 




Shift 
(%) 


TA ^^O \^ ^1 




1 9 


Q 

o 


90 
zu 




"^1 


-U.O 


\-Jr\.O^J\J 1 -LUJ 


untaggGti 


1 A 


Q 


O.iJ 


z.ou 


-u.ui 


^ 9 
-o.z 


±j!t.DOW I -LUJ 




99 


Q 


u.o 


1 '^fi 
1 .ou 


-u. yy 


9 

-u.z 


HR^ [171 


untaggGti 


9Q 

zy 


7 










1 1 V 1 lyj 


imtaggcd 


9Q 
zy 


9n 




1 ns 

1 .uo 










ou 


Q 

o 


90 
zu 




-u.oz 


Ifi 


T A ■RF' [901 


untaggGd 




9 


1 A 


07 
D.U / 


n 5i9 

-U.oZ 


-lO.O 


TAC!c;n [901 
li-lDoVJ IZUJ 


untagged 




7 





1 1 9 
J LZ 


-i.Oo 


-y .o 


TDPA 7 [91 1 




oo 


c: 
u 




n 9Q 

u.zy 






AT FPH [99 9Sl 


untagged 


Ql 9 


1 8 


Q 
o 


u.uo 


1 7^^ 


-o.o 


r»PT PWT [9*^1 


1 tagged 


Ql 9 
y i .z 


1 7 




n 91 






DPT PHT [9*^1 




Ql 9 
y ± .z 


1 7 




1 nn 






DPT PHT [9*^1 


untagged 


Ql 9 
y ± .z 


1 7 




n 9fi 






OPAT [^1 


11 4" Q fT<T ^i~l 

U. LdgjgjtJLl 


Ql 9 
y ± .z 


c: 
u 




9 14 






OPAT [Pil 


d tagged 


Ql 9 
y i .z 







9 48 






OPAT [p;l 


s tagged 


Ql 9 
y i .z 







o.Uo 






OPAT [f^l 


L, LdgjgjtJLl 


Ql 9 
y ± .z 


c; 
o 




9 fi7 

Z.U ( 






OPAL f5] 


L) LdgigtlLl 


91.2 


c; 
o 




1 ^1 






OPAT [9/11 


untagged 




1 n 

iU 




8 78 






ST.D f25l 




91.2 


29 




1.52 






SLD [25] 


c tagged 


91.2 


29 




1.51 






SLD [25] 


b tagged 


91.2 


29 




1.77 






SLD [25] 


untagged 


91.2 


29 




0.47 






DELPHI [26] 


untagged 


189 


3 




2.01 






RR AHMS fQ] 


y G [2.9, 31 


200 


7 


11,7,8(13), 
2,1(3) 


2.42 


-2.69, -1.71, -1.70, -0.49, -0.14 


-29.6, -12.0, -13.6, -1.0, -0.1 


y e [3.25,3.35] 


5 


5.33 


-3.68, -2.34, -2.51, -0.67, -0.27 


-40.4, -16.4, -20.0, -1.3, -0.3 


STAR [3] 


\y\ < 0.5 


200 


8 


11.7 


3.01 


-1.94 


-22.6 


Total 






309 




1.89 







TABLE X: As in Table HXl but for the charge-sign asymmetry Acp/p. 



Collaboration 


Properties 


(GeV) 


# 
data 


Norm. 

(%) 


2 

Xdf 


Aa- 


Shift 

(%) 


BR AHMS f9] 


y e [2.9, 31 


200 


7 


11,7,8(13), 
2,1(3) 


1.83 


-1.57, -1.00, -0.58, -0.29, 0.08 


-17.3, -7.0, -4.6, -0.6, 0.1 


y e [3.25, 3.351 


5 


3.77 


-3.10, -1.97, -2.12, -0.56, -0.23 


-34.1, -13.8, -16.9, -1.1, -0.2 


STAR [3] 


\y\ < 0.5 


200 


8 


11.7 


1.51 


-0.16 


-1.8 


Total 






20 




2.19 
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TABLE XI: As in Table |Vl but for A'g production. The CDF data are normaUzed by the unknown total cross section, which 
is therefore fitted to obtain the result a = 13.7 mb. 



Collaboration 


Properties 


(GeV) 


# 
data 


Norm. 

(%) 


2 

Xdf 


Aa- 


Shift 

(%) 


TASSO 


30] 


untagged 


14 


8 


15 


0.53 


-1.32 


-19.8 


TASSO 


31] 


untagged 


14.8 


8 




0.60 






TASSO 




untagged 


21.5 


5 




0.44 






TASSO 


30] 


untagged 


22 


5 


15 


0.79 


-0.08 


-1.3 


HRS [17] 


untagged 


29 


12 




2.61 






MARK 11 [27] 


untagged 


29 


17 


12 


0.76 


0.26 


3.1 


TPC [32] 


untagged 


29 


7 




0.41 






TASSO 


33] 


untagged 


33.3 


7 


15 


0.55 


0.26 


3.9 


TASSO 


30] 


untagged 


34 


13 


15 


0.93 


-0.05 


-0.8 


TASSO 


3^ 


untagged 


34.5 


13 




1.72 






CF,T,T,0 [341 


untagged 


35 


9 




0.46 






TASSO [31] 


untagged 


35 


13 




1.78 






TASSO [31] 


untagged 


42.6 


13 




0.88 






TOPAZ [21] 


untagged 


58 


4 




0.05 






ALEPH [28] 


untagged 


91.2 


16 


2 


0.44 


-1.52 


-3.0 


DELPHI [35] 


untagged 


91.2 


13 




0.55 






OPAL [5] 


u tagged 


91.2 


5 




0.98 






OPAL [5] 


d tagged 


91.2 


5 




0.63 






OPAL [5] 


s tagged 


91.2 


5 




0.47 






OPAL [5] 


c tagged 


91.2 


5 




1.54 






OPAL [5] 


b tagged 


91.2 


5 




1.98 






OPAL [36] 


untagged 


91.2 


16 


6 


0.46 


-0.40 


-2.4 


SLD [37] 


1 tagged 


91.2 


9 




0.71 






SLD [37] 


c tagged 


91.2 


9 




0.91 






SLD [37] 


b tagged 


91.2 


9 




1.59 






SLD [37] 


untagged 


91.2 


9 




1.17 






DELPHI [26] 


untagged 


183 


2 




11.39 






DELPHI [26] 


untagged 


189 


3 




3.81 






RR AHMS M 


y e [2.9, 3] 


200 


8 


11,7,8(13), 
2,1(3) 


1.37 


-1.24, -0.79, 0.00, -0.23, 0.07 


-13.6, -5.5, 0.0, -0.5, 0.1 


y e [3.25, 3.35] 


6 


1.81 


-0.55, -0.35, -0.91, -0.10, -0.15 


-6.1, -2.5, -7.3, -0.2, -0.2 


CDF [i4] 


\v\ < 1 


630 


48 




0.59 


-2.19 


-21.9 


STAR [3] 


\y\ < 0.5 


200 


9 


11.7 


1.06 


-1.58 


-18.5 


Total 






323 




1.03 
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TABLE XII: As in Table |Vl but for A/A production. The CDF data are normaUzed by the unknown total cross section, which 
is therefore fitted to obtain the result a = 4.5 mb. 



Collaboration 


Properties 


(GeV) 


# 
data 


Norm. 

(%) 


2 

Xdf 




Shift 

(%) 


TASSO [30] 


untagged 


14 


3 


20 


0.17 


0.05 


1.1 


TASSO [30] 


untagged 


22 


4 


20 


0.43 


-0.60 


-12.0 


HRS [38] 


untagged 


29 


12 




0.93 






MARK II [39] 


untagged 


29 


15 




0.71 






TASSO [33] 


untagged 


33.3 


6 


15 


0.73 


-1.36 


-20.5 


TASSO [30] 


untagged 


34 


6 


20 


0.48 


-0.96 


-19.2 


TASSO [40] 


untagged 


34.8 


9 


9 


1.88 


0.26 


2.3 


CELLO [34] 


untagged 


35 


7 




1.02 






TASSO [20] 


untagged 


42.1 


4 


9 


0.22 


-0.27 


-2.4 


ALEPH [28] 


untagged 


91.2 


16 


4 


0.34 


-1.35 


-5.4 


DELPHI [41] 


untagged 


91.2 


7 




1.33 






OPAL [5] 


u tagged 


91.2 


5 




0.22 






OPAL [5] 


d tagged 


91.2 


5 




0.20 






OPAL [5] 




91.2 


5 




0.43 






OPAL [5] 


c tagged 


91.2 


5 




1.20 






OPAL [5] 


b tagged 


91.2 


5 




0.32 






OPAL [42] 


untagged 


91.2 


12 




2.61 






SLD [37] 


1 tagged 


91.2 


4 




3.40 






SLD [37] 


c tagged 


91.2 


4 




1.50 






SLD [37] 


b tagged 


91.2 


4 




2.70 






SLD [37] 


untagged 


91.2 


9 




0.67 






DELPHI [26] 


untagged 


183 


3 




6.42 






DELPHI [26] 


untagged 


189 


3 




6.11 






CDF [14] 


k;l < 1 


630 


34 


10 


1.27 


-3.18 


-31.8 


STAR [13] 


\y\ <o.5 


200 


9 


11.7 


5.84 


-6.81 


-79.5 


Total 






188 




1.45 
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^ n+X. V5=200 GeV PP (PP for CDF) ^ +X, Vi=200 (630 for CDF) GeV 




pp p/p +X, V.v=200 GeV 




P^(GeV) 



FIG. 1: Comparison of the calculation with measurements of the invariant differential cross sections for inclusive production 
in PP collisions at ^/s = 200 (630 for CDF) GeV. Note that Ed^a/dp^ is averaged over rapidity using the operation (l/(j/max — 
2/min)) jy^°'^ dy. For 7r*, the curves for which —0.35 < y < 0.35 and —0.5 < y < 0.5 overlap, as do those for which 3.25 < y < 
3.35 and y = 3.3. 
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pp %-+X, BRAHMS (2.9<.v<3), V.v=200 GeV 




pp^IC +X, BRAHMS (2.9<y<3), V.s=200 GeV 




pp p/p +X, BRAHMS (2.9<y<3), V.s=200 GeV 
lo S j 




FIG. 2: As in Fig. [T]for the BRAHMS data for which 2.9 < j/ < 3 (labeled "Default"). Also shown is the calculation when 
mass effects are neglected (the dotted curve labeled "rrih = 0") and when the ratio k — Mj/s is increased to 4 (lower solid 
curve) and decreased to 1/4 (upper solid curve). In the case of vr*, the = curve cannot be seen because it overlaps with 
the "Default" curve. 




FIG. 3: As in Fig. but for the STAR data. 
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PP ■ 



■ n'+X, BRAHMS (2.9<.v<3), V.v=200 GeV 



Ed^a 



dp 

[mb(GeV)"^] 



10 

lo' 

lo" 

10-' 
10-^ 

10"' 

10-^ 

10° 

10-' 
10-^ 



1 1 1 1 1 1 1 1 1 1 1 1 


II 


- 




— AKKII 






AKK 






DSS 






- HKNS 










" , , , 1 , , , 1 , , , 







12 3 4 

pp^IC +X, BRAHMS (2.9<y<3), V.s=200 GeV 




1 2 3 4 5 

PP p/p +X, BRAHMS (2.9<y<3), V.s=200 GeV 




2 3 

Pr(GeV) 

FIG. 4: As for the "Default" cases in Fig. [2] (now labeled "AKKQ8"), including also the calculations using the AKK jjjiii], DSS 
[m and HKNS [l3| FF sets with the hadron mass set to zero. 



pp n-+X, STAR (-0.5<3;<0.5), Vj=200 GeV PP ^ K° +X, STAR (-0.5<>'<0.5), V.?=200 GeV 




pp p/p +X, STAR (-0.5<y<0.5), Vi=200 GeV ^ A/A +Z f-0.5<y<0.5), V.?=200 GeV 




FIG. 5: As in Fig. H] but for the STAR data. Also shown, for A/A, is the calculation using the DSV ^] FF set. 



pp ^ A n~+X. Vi=200 GeV 



Ed^a 



[mb(GeV)"^] 10 ' 




STAR (-0.5 <>'<0.5) 
BRAHMS (2.9<>'<3) 
BRAHMS (3.25<v<3.35) 



I — STAR (-0.5<v<0.5) 
BRAHMS (2.9<v<3) 
BRAHMS (3.25<v<3.35) 



9 1( 



pp ^ A IC+X, Vs=200 GeV 




12 3 4 

pp A j)/p+X, Vi-200 GeV 




3 4 

P^(GeV) 



FIG. 6: As in Fig. [T] but for the charge-sign assymetry Ac/i^ 
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pp A 7t"+X, BRAHMS (2.9<v<3), Vi=200 GeV 




pp ^ A tr+X, BRAHMS (2.9<y<3), Vi=200 GeV 




pp ^ A jy/p+X, BRAHMS (2.9<y<3), V.s=200 GeV 



10 
10^ 
10-' 

10-= 

■A 

10 
10-" 

10-^ 
10-^ 
10-^ 
10-' 




P^(GeV) 



FIG. 7: As in Fig. |6]for the BRAHMS data for which 2.9 < y < 3 (labeled "Default"). Also shown is the calculation when 
mass effects are neglected (the dotted curve labeled "mji = 0") and when the ratio k — Alj/s is increased to 4 (lower solid 
curve) and decreased to 1/4 (upper solid curve). In the case of Acvr*, the ruh = curve cannot be seen because it overlaps 
with the "Default" curve. 
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pp A^7i"+X, STAR (-0.5<y<0.5), V.s=200 GeV 
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FIG. 8: As in Fig. □ but for STAR kinematics. 
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pp A%-+X, BRAHMS (2.9<v<3), Vi=200 GeV 
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FIG. 9: As in Fig.H but for Acft*. 
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pp A^7i"+X, STAR (-0.5<y<0.5), V.s=200 GeV 
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pp^A jy/p+X, STAR (-0.5<v<0.5), Vi=200 GeV 
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FIG. 10: As in Fig. H but for STAR kinematics. 
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pp^fC +X (2.9<y<3), M=p^, ^s=200 GeV 
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FIG. 11: Contributions to the production from fragmentation of the initial state protons' valence quarks and sea partons, 
for BRAHMS kinematics with 2.9 < j/ < 3. In the notation of Eq. (|62|l . "u — m" refers to the contribution from u valence 
fragmentation FuD'^ — FuD'i etc. The ratios are stacked, i.e. for a given pr value, the distance on the y-axis from zero to 
the first curve gives the u — u contribution, from the first to the second curve the d — d contribution, and from the second curve 
to 1 the sea parton contribution. 



pp it+X ('-0.5<y<0.5), M=p^, V.s=200 GeV 

1| ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' 

8- 

g+2u+2d+2s+... 
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FIG. 12: As in Fig. [HI but for STAR kinematics. 
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pp ^^^~+X f2.9<.v<3), Vi=200 GeV 

0.002 1 — . — , — , — , — ^ — . — , — , — , — ^ — , — . — , — r 
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2 2.5 3 3.5 



FIG. 13: Contributions to the charge-sign assymetry Acft for BRAHMS kinematics for which 2.9 < j/ < 3 from fragmentation 
of the initial state protons' valence quarks. In the notation of Eq. (|62[) . "it + it" means FuD^"^ + FuD^"^ etc. The cross 
section is given by the difference between the two curves. 
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pp An~+X f-0.5<j'<0.5), Vs=200 GeV 
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FIG. 14: As in Fig. [H but for STAR kinematics. 
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e e -^Tt +X 



e e ^K^+X 



1 do 10 
a dx jQ-i 



TASSO (V.v=14GeV) 
TASSO (Vi=34 GeV) 
OPAL(Vi=91.2GeV) 
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FIG. 15: Comparison of the calculation with measurements of the invariant differential cross sections for inclusive /i^ production 
in e^e~ collisions. For clarity, the TASSO data at ■\/s = 34 GeV and the OPAL data at ^/s = 91.2 GeV have been multiplied 
by 0.1 and 0.01 respectively, as well as their corresponding calculations. 




FIG. 16: As in Fig. [15] for the TASSO data at = 
effects are neglected (the dotted curve labeled "mj, — 
and decreased to 1/4 (upper solid curve). 



14 GeV (labeled "Default"). Also shown is the calculation when mass 
0") and when the ratio k = M'j / s is increased to 4 (lower solid curve) 
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e*e n+X, OPAL, a/.s=91 GeV 




e*e p/p+X, OPAL, Vi=91 GeV 




X 

FIG. 17: As in Fig. [16] but for the OPAL data at ^ = 91.2 GeV. 
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i~+X, Vi-=91.2GeV 



e*e K° +X, V.s=91.2 GeV 
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FIG. 18: As in Fig. 1151 but for all fitted -^s — 91 GeV data and witli the corresponding calculations from other FF sets. 
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e*e n+X, TASSO, V.s=14 GeV 



e*e K° +X, TASSO, V.s=14 GeV 
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FIG. 19: As in Fig. 1181 but for the TASSO data at -^s — 14 GeV data and with the corresponding calculations from other FF 
sets. 
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FIG. 21: OPAL quark tagging probabilities for tt*. 
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FIG. 22: OPAL quark tagging probabilities for K 
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e e —>uu—> 



p/p+X, OPAL, Vi=91 GeV 



e^e -^cc^ p/p+X, OPAL, Vi=91 GeV 
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e*e -^ss^ p/p+X, OPAL, ^s=9l GeV 



10" 




FIG. 23: OPAL quaxk tagging probabilities for p/p. 
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ee ->dd^Kl +X, OPAL, Vi=91 GeV 



e^e -^cc^K^ +X, OPAL, ^5=91 GeV 
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FIG. 24: OPAL quark tagging probabilities for Kg. 
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FIG. 25: OPAL quaxk tagging probabilities for A/A. 
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H=TC, i=u, M=912 GeV 
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FIG. 26: The FFs for tt* at M/ = 91.2 GeV. 
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H=Kr, i=u, M=91.2 GeV 



H=K-,i=c, M =91.2 GsV 
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FIG. 27: The FFs for at M/ = 91.2 GeV. 
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H=p/p, i=u, M=91.2 GeV 



H=p/p, i=c, M=9\.l GeV 
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FIG. 28: The FFs for p/p at My = 91.2 GeV. 
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FIG. 29: The FFs for K°s at Mj = 91.2 GeV. 
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FIG. 30: The FFs for A/A at M/ = 91.2 GeV. 



56 



H=A/A, i=c, M=9\.2 GeV 
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FIG. 31: The quark FFs for Ac7r=^ at M/ = 91.2 GeV. 



H=A IT, i=M, M =91.2 GeV 




FIG. 32: The u quark FF for Aei^=^ at M/ = 91.2 GeV. 



